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The temporal evolution of step-edge fluctuations under electromigration conditions is analysed 
using a continuum Langevin model. If the electromigration driving force acts in the step up/down 
direction, and step-edge diffusion is the dominant mass-transport mechanism, we find that significant 
deviations from the usual t 1 / 4 scaling of the terrace-width correlation function occurs for a critical 
time t which is dependent upon the three energy scales in the problem: ksT, the step stiffness, 
7, and the bias associated with adatom hopping under the influence of an electromigration force, 
±AJ7. For (t < t), the correlation function evolves as a superposition of t 1 ^ 4 and i 3//4 power laws. 
For t > r a closed form expression can be derived. This behavior is confirmed by a Monte-Carlo 
simulation using a discrete model of the step dynamics. It is proposed that the magnitude of the 
electromigration force acting upon an atom at a step-edge can by estimated by a careful analysis of 
the statistical properties of step-edge fluctuations on the appropriate time-scale. 

PACS numbers: 05.40.-a, 66.30.Qa, 68.35.Ja 



I. INTRODUCTION 

During the past decade, continuum models and dis- 
crete lattice simulations have been applied to understand 
direct imaging observations of the thermal fluctuations of 
step edges in which the step position is monitored as a 
function of time [H, 0| ■ Of particular interest has been 
the study of the dynamics of the equilibration of terrace 
width distributions where the temporal evolution of step 
edge fluctuations are driven by the exchange of atoms 
between the step and the adjacent terrace and/or bymo- 
tion of adatoms along the step edge itself [H 0, H [1,0, [H • 
It is well known that the position of the step edge, as de- 
scribed by its temporal correlation function, has a time 
dependence that scales as a power law with an expo- 
nent characteristic of specific atomic processes driving 
the step fluctuations; t@ . In cases where mass transport 
at the step is dominated by diffusion of atoms along the 
step edge; (3 — 1/4. When mass transport proceeds via 
exchange of atoms between the step edge and the terrace 
(3 = 1/2. Careful experiments allow a crossover from i 1 / 4 
to i 1 / 2 scaling to be observed as a function of the sam- 
ple temperature @ . Further, experimental measurement 
of the correlation functions have been used to determine 
thermodynamic properties of the steps, such as the step 
stiffness 111 0]. 

In this paper we investigate how the scaling of the step 
edge fluctuations is changed by the presence of an electro- 
migration force [l(| acting upon atoms diffusing along the 
step edges. The primary motivation for this study is the 
possibility of using measurements of these changes to ob- 
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tain information about the electromigration force itself. 
In conducting materials, a surface electromigration force 
can be generated by passing an electrical current through 

the sample m E3, Hi 0, til M Gl M M M MM- 

In terms of a simple discrete model, the presence of the 
electromigration force introduces a small bias in the dif- 
fusion of atoms at the step edge, parallel to the current 
(and field). By convention, this bias can be expressed as 
an energy difference between atoms diffusing parallel or 
anti-parallel to the current AU = Z*eEa± where E is 
the electric field applied to the sample, Z* is the effec- 
tive valence of adatoms and a± is the lattice spacing per- 
pendicular to the step edge. A characteristic property of 
surface electromigration is that the electromigration bias, 
AU, is several orders of magnitude smaller than the other 
energy scales in the problem: 7a j_, where 7 is the step 
stiffness and the energy associated with thermal fluctua- 
tions; fcflT. This suggests that thermal fluctuations will 
completely dominate the short-time behavior of the step 
fluctuations with the effect of the electromigration bias 
emerging only on much longer time scales. Nevertheless, 
such time scales (of the order of seconds) are accessible to 
experiment offering the possibility that the observation of 
step fluctuations under electromigration conditions may 
allow us to obtain quantitative information about the 
magnitude of the force itself; a quantity that, to date, 
has been hidden from experimental study. 

This paper is organized as following. In section 2 we 
present a continuum theory of step edge fluctuations un- 
der electromigration conditions. In section 3 we test the 
fidelity of the theory by showing the results of a Monte- 
Carlo simulation of the temporal evolution step correla- 
tion function. Concluding remarks are contained in sec- 
tion 4. 
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II. THEORY 

In order to describe the dynamics of a step-edge evolv- 
ing under electromigration conditions we employ the 
usual Langevin formalism where each degree of freedom 
diffuses towards lower energy with a velocity that is pro- 
portional to the energy gradient, subject to random ther- 
mal fluctuations. The position of step edge is described 
by it's edge profile y(x,t) where the x-axis is oriented 
parallel to the step edge. y(x,t) is the position of the 
step edge at x and at time t. 

In this paper we consider the limiting case where the 
step motion occurs most easily by adatom diffusion along 
the step edge itself. Adatom exchange between the step- 
edge and the adjacent terraces, via attachment and de- 
tachment, is neglected. It is well known that atomic diffu- 
sion along a step edge is driven by the step curvature [23| , 
which generates a flux 



kT 



(1) 



where J c is the curvature-driven flux, D s = ctfJrh is 
the surface diffusion constant, an is the atomic diameter, 
S s = aj_ the width of an atom perpendicular to the step, 
Th is the average time between hopping events and 7 is 
the step stiffness. V s k is the gradient of the curvature 
along the step edge. Mass conservation determines the 
normal velocity of the step edge, 



dy 
at 1 



-fiV s • J c 



(2) 



where Q is the area of a single atom and h is normal to 
the step edge. 

The inclusion of a thermal noise term, rj and lineariza- 
tion of the above equation leads to the well-known equa- 
tion of motion for a step edge [|| 



d T hl d 4 



dt kT dx 4 



where we have defined 



y(x,t) = J](x,t). 
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We model the effects of electromigration by adding to 
this equation of motion a term generated by a constant 
force, F = Z*e~E, felt by atoms at the step edge, which 
arises from the application of a electric field, E, to the 
material oriented parallel to the y-axis. Z* is the effective 
valence of an atom at the step edge. The electromigration 
force generates an additional flux, 
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where E\ \ is component of the electric field along the step 
edge. y x indicates an x derivative. 



The energy of any step configuration, relative to the 
energy of the straight step, is determined by the step stiff- 
ness 7 and the electromigration force, F, felt by atoms 
at the step edge. If the force acts perpendicular to the 
step edge (i.e. when F > 0, the force acts in the +y, 
or step down, direction) then, for small fluctuations, one 
can linearize the stochastic equation of motion for the 
step edge to obtain the following Langevin equation for 
the step dynamics, 



dt kT dx 4 kTa\\a± dx 2 



y(x,t) = Tj(x,t) (6) 



where an and a± are the lattice parameter parallel and 
perpendicular to the step edge. The noise term, 77, de- 
scribes the thermal fluctuations and is correlated be- 
cause, in our model, of the random hopping of atoms 
occurs only between between adjacent step edge sites. 



{r)(x,t)r)(x' ,t')) 



2T h S(t-t')^5(x-x') (7) 



To first order, the electromigration force does not al- 
ter the correlation properties of the noise so that in the 
single-jump regime equation [4] remians valid. 

Before proceeding, for notational simplicity we rewrite 
eqn. \E\ as 



0_ 

dt 



a dx~ 4Taq ^ 



y{x,t)=r)(x,t) (8) 



where 



a 
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-F=i/ m < 9 > 



The parameter q c depends only on the magnitude of the 
electromigration force and in eqn. [6] the =F denotes the 
force acting in the step down(-) or step up (+) direction. 

In the case where there is no random thermal noise 
(i.e. T — > 0), eqn. [6] predicts that the amplitude of a 
small sinusoidal fluctuation of the step-edge profile, with 
wavevector q, evolves according to the following disper- 
sion relation, 



aq' 



(10) 



This is the same dispersion relation as the one that de- 
scribes step flow under growth conditions [13, HE Hfi| • If 
the electromigration force acts in the step up (+) direc- 
tion the step fluctuations are linearly stable. For a force 
acting in the step-down direction there exists a linear in- 
stability which initiates the well-known phenomenon of 
electromigration-induced step- wanderin g fo r fluctuations 
with wavevector smaller than q c [27], I28l. l29l| . The critical 
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wavevector is an important parameter in our Langevin 
model that, from eqn. [9] is determined by the ratio of 
the electromigration force and the step stiffness. 

In order to determine the statistical properties of the 
solutions of eqn. [6] we take the usual approach and first 
derive the Green's function for the problem using stan- 
dard Fourier transform methods, 



/+oo 
e - a{k * +q lk*){t-t') e ik(x- X >) dk 
-oo 

(11) 

In terms of the Green's function, the displacement of the 
step at time t is: 



/+oo ft ' 
/ G{x,t"\x / ,t)r 1 (x , ,t , )dx'dt 
-oo Jt 

We can now compute the time correlation function of the 
step edge, G(t), after time t = t" — t' has elapsed, 



/(*)=< (»(*,*")-!/(*. *0r > 



(13) 



Substituting eqn. QT] into [12] and using the correlation 
properties of the noise (eqn. [7]) we obtain we find, after 
some calculation, 



9±(t) 



<^ i '/.") 



or, using substitution, 



(14) 



1/4 
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1/4 (kTyr h \ ■ 2i 



(15) 



~ 2u e^ ^Adu 



o (u 2 ± ^/ak£) 



When there is no electromigration force (i.e. \F\ = 
0, q c = 0), the integral in eqn. \W\ is clearly time- 
independent and we recover the result derived by Bartelt 
et al. [3J where the g(t) the step edge fluctuations evolves 
with the well-known t 1 / 4 scaling law characteristic of 
step-edge limited diffusion. 



Bm ftt (t) = ft ,(t)=tV* (i^DllH 



(16) 

It is helpful to re-express eqn. [TC] in terms of the average 
time, To, that it takes for the mean-squared width of an 
initially straight step to reach a value equal to 5o( r o) = 
a\ (i.e. one square lattice spacing): 



9o(t) 
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where 



TO = Th 



\2kTa ) \ 27rcV4) 



1<A 



(18) 



When the electromigration force is finite ( i.e. \F\ > 0), 
it is apparent from eqn. [16] that this scaling behavior is 
modified by the appearance of an explicit time depen- 
dence in the integral. This can be seen more clearly by 
defining a critical time, t, and a rescaled time, C = t/ T i 
where, 



aq A c 



\Fa± J \2fcTa|| 
Then, eqn. [16] can be rewritten as 
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9±{t) 
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where 

i±(0 



(20) 



(21) 



2V4r(3/4) 
+°° 1 / 
(u 2 ±Vt) { 



I± is a universal function of the rescaled time, £ and 
is normalized such that I±(t — > 0) = 1. The integral 
appearing in eqn. [22] is easily evaluated numerically and 
is shown in fig. 1 in which we display I± plotted as a 
function of the rescaled time £ (i.e. time is expressed in 
units of t) . The solid curves show I± obtained for F > 
(step up) , F < (step down) electromigration forces 
(ForF = 0,/±(C)=l). 

From fig. 1. it is apparent that g±(t) (eqn. l2"0"]) deviates 
very significantly from the i 1 / 4 scaling behavior of go(t) 
feqn. [T7|) as t approaches r (C — > 1). These deviations be- 
gin to appear at earlier times, t ~ 0.1t, when the effect of 
the force on the evolution of the step fluctuations begins 
to be felt. This can be seen more clearly in fig. 2 which 
displays the correlation function of a step plotted as a 
function of time for To = 5s and r = 10 +4 s. The values 
for these parameters were chosen so that t/tq ~ 10 4 , a 
ratio typical of accelerated electromigration experiments. 
As noted above, deviations from i 1 / 4 scaling start to ap- 
pear when t ~ 0.1t = 100s. The results shown in figs. 1 
and 2 have a simple qualitative interpretation; the short 
time behavior of the step fluctuations (t -C t) is com- 
pletely dominated by the thermal fluctuations and the 
effect of the electromigration bias emerges only at later 
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times. Such behavior is typical of the dynamics of diffu- 
sion driven by weak external forces. 

It is instructive to perform a power series expansion of 
the integral about ( = such that eqn. [20l becomes. 
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The expansion coefficents can be obtained analytically, 



a i 



er(§)2 



0.3487, 



20 = °- 150 °' 



0.2 0.4 0.6 0.8 1 



(23) 

Shown as the dashed lines in fig. 1 are the results of this 
series approximation for I±{Q (eqn. [52]), evaluated up 
to, and including the terms linear in time. Clearly, this 
truncated expression is a reasonable approximation for 
t < 0.4t. 



FIG. 1: The integral function (eqn. [22)) plotted as a 

function of the rescaled time £ = t/r for no electromigration 
force (F — 0), the electromigration force in the step down 
direction (F > 0) and in the step-up direction (F < 0). 
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FIG. 2: The time corrleation function, g(t), of the step- 
edge distribution predicted by the continuum Langevin model 
(eqns. |2"01 and |2"2")) . plotted as a function of time for to = 5s 
and r = 10 4 s. 



III. SIMULATION 

In order to test the predictions of the continuum 
Langevin model described above, we developed a Monte 
Carlo simulation of step edge fluctuations in the pres- 
ence of an external force. In this model, atomic diffu- 
sion was restricted to the step edges with atoms jump- 
ing between adjacent step sites only. Only nearest- 
neighbor interactions on a square lattice were permit- 
ted (dj_ = an = a = 1) and were modelled by a single 
bond energy e. In order to describe the electromigration 
force, the atomic jump probabilities for motion parallel 
and anti-parallel to the force were biased by a small en- 
ergy differential AU. In terms of the electromigration 
force, F, and the lattice spacing perpendicular to the 
step edge, a; AU — ±Fa. 

Simulations were performed for steps of length I — 
10000a|| fluctuation on a square lattice. Periodic bound- 
ary conditions were employed. The bond energy was set 
to £ = 2.0/cbT and the magnitude of the electromigra- 
tion bias was \AU\ — O.OlfcsT, a factor of 10 3 smaller 
than the typical binding energy of an atom to the step 
edge. This value was chosen to generate statistically sig- 
nificant deviations from the (no-force) i 1 / 4 scaling within 
reasonable simulation times. If e = O.leV and it is as- 
sumed that a ~ 1.5 A (typical of metals) then this bias 
corresponds to an electric field with a magnitude of or- 
der lOOOI^cm -1 acting on an atom with effective valence 
of Z* ~ ±10e. In actual accelerated electromigration 
experiments, a field of 0.1 — WcmT 1 is typical. 

Figure 3 shows the results of the simulation where the 
correlation function of an isolated step is plotted as a 
function of the time measured in Monte-Carlo steps per 
step-edge site (MCS). Shown is g(t) obtained when the 
electromigration force acts in the step-up and step-down 
directions and when AU — (i.e. no electromigration 
force is acting at the step edge). We define one Monte 
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FIG. 3: Results of the a kinetic Monte Carlo simulation where 
the correlation function of an isolated step is plotted as a 
function of the time (measured in Monte-Carlo steps per step- 
edge site, the lattice spacing is a = 1). Shown is g(t) obtained 
when the electromigration force acts in the step-up and step- 
down directions and when AU = (i.e. no electromigration 
force is acting at the step edge). The curves shown were 
obtained by averaging the data from 200 randomly generated 
replicas after each was equilibrated for 10 5 Monte Carlo steps. 



Carlo step to be equal to the average time needed for ev- 
ery atom at the step edge to attempt a hop. The results 
shown in fig. 3 were obtained by averaging the data from 
100 randomly generated replicas after each was equili- 
brated for 10 5 Monte Carlo iterations per site. Compar- 
ing the simulation results (fig. 3.) to the predictions of 
the Langevin theory (fig 4.) it is apparent that the qual- 
itative features of the continuum theory are reproduced 
by the Monte Carlo simulation. These same data are pre- 
sented in the form of a log-log plot in fig. 4. In the ab- 
sence of the force (F = , AU = 0), a least-squares fit of 
the no- force simulation data shows that go(t) is very well 
fit by a ^ power-law where [3 = 0.25 ± 0.01. Therefore, 
when there is no electromigration force present in the 
simulation, the correlation function of the step evolves 
according well-known t 1 / 4 scaling, as predicted by the 
Langevin analysis (eqn. [T7|) . By least squares fitting the 
simulation results to eqn. [T7] we obtain a value of tq = 5 
MCS. 

We now consider the results of the simulation obtained 
for a finite electromigration force, also shown in figs. 3 
and 4. Equation [2U suggests that the value of r can be 



FIG. 4: A log-log plot of the simulation data shown in fig. 3. 
The dashed line shows the best fit of a power law, (t/ro) 13 , to 
the no-force (F = 0) data; (3 = 0.25, r = 5 MCS. 

extracted from the simulation results by considering the 
scaling of the difference between the correlation functions 
for the force in the up and down step directions: 



A -2 
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~t 




+ o 




T 




T 



(24) 



For the simulation results, this normalized difference is 
plotted in fig. 5. The behavior of this quantity is well 
fit by the leading order term in eqn. [24] from which we 
obtain a value of r = 62000 ± 10000 MCS. 

For comparison with the fits to the continuum 
Langevin model , we can estimate r from the microscopic 
parameters employed in the discrete Monte Carlo model 
used to generate the simulation data. Combining eqn. [TO] 
with the step stiffness appropriate for our model 

^"j = ea ± (25) 
2/cTa|| 2k B Ta\\ v ' 

we obtain an estimate of r in units of MCS: 

T = 2 ( T ±]( k -f-\ 2 ( ^ ) (26) 
\tJ \AU J \2k B TaJ 1 ' 

The ratio of the hopping time in the Langevin model , 
Tfi to the time between hopping attempts in the Monte 
Carlo simulation r a can be obtained by monitoring the 
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FIG. 5: For the simulation results shown in figs. 3 and 4; a 
log- log plot of the normalized difference (eqn. I24p plotted as a 
function of time (MCS). The behavior of this quantity is well 
fit by the leading order term in eqn. [24] (t 1 ^ 2 , dashed curve) 
from which we obtain a value of r = 62000 ± 10000 MCS. 



success rate of hops between adjacent lattice sites in the 
simulation. We find that r a /r^ = 3.6 ±0.1. Thus our 
estimate for the value of t in the Monte Carlo simu- 
lations, used to generate the data shown in figs 3-5, is 
r = 62000 ± 10000 MCS. Clearly, the agreement between 
the continuum Langevin theory (r = 71000 MCS ) and 
the microscopic model (r = 62000 ± 10000 MCS) is rea- 
sonable. Finally we note that in the high temperature 
limit (fcsT 3> e) the ratio of the electromigration bias to 
the binding energy at a step edge, e is related directly to 
t that would be obtained from experiment: 



where Th can be determined from eqn. 1181 if the step 
stiffness is known. 



IV. CONCLUSIONS 

The temporal evolution of step-edge fluctuations un- 
der electromigration conditions has been analyzed using 
a continuum Langevin model for the case where diffu- 
sion is limited by mass transport along the step edge. 
We find that the presence of the electromigration force, 
felt by atoms at the step edge, causes deviations from 
the usual t 1 ! 4 scaling of the terrace-width distribution 
driven by thermal fluctuations alone. We have identified 
a critical time r that is a function of the three energy 
scales in the problem: fc^T, the step stiffness, 7, and the 
bias associated with adatom hopping under the influence 
of an electromigration force, ±A{7. For (t < r), the cor- 
relation function evolves, to a good approximation, as a 
superposition of t 1 / 4 and t 3 / 4 power laws. For all r a 
closed form expression was derived. This behavior was 
confirmed by a Monte-Carlo simulation using a discrete 
model of the step dynamics. Finally we propose that the 
magnitude of the electromigration force acting upon an 
atom at a step-edge could be determined directly by care- 
ful measurement and analysis of the statistical properties 
of step-edge fluctuations on the appropriate time-scale. 
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